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Abstract. We consider bounds on codes in spherical caps and related prob- 
lems in geometry and coding theory. An extension of the Delsarte method is 
presented that relates upper bounds on the size of spherical codes to upper 
bounds on codes in caps. Several new upper bounds on codes in caps are de- 
rived. Applications of these bounds to estimates of the kissing numbers and 
one-sided kissing numbers are considered. 

It is proved that the maximum size of codes in spherical caps for large 
dimensions is determined by the maximum size of spherical codes, so these 
problems are asymptotically equivalent. 

1. Introduction 

The subject of this paper is codes in spherical caps, i.e., packings of a metric ball 
on the surface of the sphere in W 1 (a spherical cap) with metric balls (caps) of a 
smaller radius. Codes in spherical caps are related to more familiar spherical codes 
and find a number of interesting applications in both classical and recent works. 

Spherical cap codes have been used to derive an asymptotic upper bound on 
the maximum size of spherical codes and a bound on the packing density A„ of 
the n-space by equal spheres, see Sidelnikov [24] . Kabatiansky and Levenshtein [TU] 
and Levenshtein [12] . More recently they have been used to derive upper bounds 
on the size of binary constant weight codes, see Agrell, Vardy and Zeger pQ. Even 
more recently, estimates of the maximum size of codes in a spherical cap have been 
used together with an extension of Delsarte's method to derive new estimates of the 
kissing number k(n) in small dimensions. In particular, a long-standing conjecture 
that fc(4) = 24 was solved in [THl E] and a related problem for "one-sided kissing 
numbers" was solved in dimension 4 in [15] . 

In this paper we focus on the study of spherical cap codes rather than on their 
applications in related geometric problems of coding theory. More specifically, we 
study bounds on the size of spherical cap codes with a given angular separation and 
their relation to spherical codes. 
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In Section [3] we recall a few known bounds on the size of codes in a spherical 
cap and spherical strip. In Section U we formulate an extension of Delsarte's bound 
on the size of spherical codes to cover the case of spherical caps. As usual, the 
polynomial involved in the computation of the bound must be expandable in a 
linear combination of Gegenbauer polynomials with nonnegative coefficients. A 
new condition in the theorem relates the values of the polynomial to constructions 
of codes in a spherical cap. The method described was used in [14j to prove that 
k(A) = 24. This link serves as an additional motivation for studying spherical cap 
codes. 

In Section[5]we show that in the case of large code distance, the size of a code in a 
spherical cap can be exactly expressed via the size of codes on the entire sphere. The 
result is used to relate the size of spherical cap codes with the kissing number k(n). 
We consider examples of small dimensions n = 3, 4 and illustrate the application 
of the extended Delsarte's method to the derivation of the values of k(n) in these 
cases. 

In Section [6] we derive a new bound on the size of spherical cap codes that relies 
on a transformation from codes in caps to codes on the hemisphere. In the same 
section we also address the problem of the maximum size of spherical cap codes in 
the case of large dimensions. A common perception in coding theory, originating 
with the asymptotic results of [10] is that codes in spherical caps are analogous to 
constant weight binary codes (i.e., codes formed of vectors with a fixed number of 
ones). Constant weight codes possess a rich combinatorial structure related to the 
properties of the Johnson graph [7]; however, no similar theory has arisen for the 
spherical case. We provide an explanation of this by showing that the asymptotic 
problem of constructing spherical cap codes is equivalent to the analogous problem 
for codes on the entire sphere. 

Section [7] is devoted to a particular case of cap codes, namely codes in hemi- 
spheres. We derive an upper bound on the size of such codes and use it to derive 
estimates of a parameter closely related to k(n), the so-called one-sided kissing 
number B(n). We derive estimates of B(n) for n = 5,6, 7,8 and make a conjecture 
about the exact values for some of these cases. 

In Section [8] we use the method of Section [7] to derive another upper bound 
on spherical cap codes that often improves the result of Section [6] We show in 
examples that the bounds derived in this paper are sometimes better than the 
previously known results. In Section [5] we present a general approach to bounding 
the size of codes in spherical caps that combines several features of the methods 
introduced earlier in the paper. We conclude with a brief discussion of applications 
of the bounds on cap codes. 

2. Notation and preliminaries 

Let iS" -1 be a unit sphere in n dimensions and let e„ = (0, ...,0,1) be the 
"North pole." Let < -0 < < 90°. Denote by 

Z(n, [ip,4>]) = {x G S" 1_1 | cos0 < (x,e n ) < cos-0} 

a strip cut on the sphere by two planes perpendicular to the vector e n . In particular, 
Cap(n, 4>) — Z{n, [0,^]) is a spherical cap with angular radius <f> drawn about e n . 
A finite subset C C Z is called a code. Below by dist(-, •) we denote the angular 
distance between two points on the sphere. If a code C has minimum angular 
separation 8, i.e., satisfies dist(a;i, £2) > for any two distinct points x±,X2 G C, 
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we call it a 9- code. Let 



A(n,6, hM]) 



max I CI 

CC2(n,[i*]) 



C a 0-codc 



be the maximum size of a 9-code in the strip Z. For spherical caps we will write 
A(n, 9, 4>) instead of A(n, 9, [0, cf>]) and use a separate notation B(n, 9) :— A(n, 9, 7r /2) 
for codes in the hemisphere S+ :— Cap(n, */2). In the case of cj) = ir (the entire 
sphere) we will call C a spherical 9-code and use the notation A(n, 9) to denote its 
maximum possible size. 

The quantity k(n) = A(n, ^/z) is equal to the number of nonoverlapping unit 
spheres that can touch the sphere S*"" 1 and is called the kissing number \a dimension 
n. The problem of finding or bounding k(n) has a rich history [4l [Til [T9l 120] . 

Likewise, the quantity B(n) = B(n,*/3) is called the one-sided kissing number. 
The one-sided kissing number problem was considered recently in [US]. B(n) has 
the following geometric meaning. Let H be a closed half-space of the n-dimensional 
Euclidean space. Suppose that S is a unit sphere in H that touches the supporting 
hyperplane of H. The one-sided kissing number B{n) is the maximal number of 
unit nonoverlapping spheres in H that can touch S. 

The function A(n, 9) has received considerable attention in the literature. There- 
fore, one possible avenue of studying spherical cap codes is to map them on the 
sphere or hemisphere and relate them to spherical codes. In this paper, we rely 
on a number of mappings between spheres, spherical caps, and spherical strips to 
estimate the maximum size of a code in a spherical cap. Some of them have been 
used earlier in the literature while the others have not been emphasized in the con- 
text of estimating the code size. The main problem addressed here is to design the 
mappings so that the distance between the images of two points in the domain can 
be bounded in terms of the original distance. One often-used map is the orthogonal 
projection LT n which sends the point x £ S n ~ l along its meridian to the equator 
of the sphere, i.e., the set of points on S 71-1 with x n = (n n is defined on S 71 ' 1 
without the North and South poles). Below we use the notation S n ~ 2 to refer to 
the equator of the sphere S 71 ^ 1 . 

Throughout this paper we use the function ui(9, a, (3) defined by 



In the case of a = (3 we write us {9, a) instead of u>(9, a, a). This function describes 
the change of the distance between two points on S 71 ^ 1 which are a and (3 away 
from e n and 9 away from each other under the action of II„ . 



Several estimates on the size of spherical cap codes have previously appeared in 
the literature. They connect the maximum size of codes in a spherical cap, and 
more generally, in a spherical strip and on the entire sphere. 

3-A. Spherical cap codes and spherical codes. Let m(n, d) be the maximum 
number of points in a unit ball in W l that lie at Euclidean distance d or more apart. 
Bounds on A(n, 9, <fi) are given in the following theorems. 

Theorem 1. (Sidelnikov [24], Levenshtein [IT] ) 



cosiv(9, a, (3) 



cos 9 — cos a cos (3 



sin a sin (3 



3. Spherical strip (cap) codes and spherical codes 
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The proof is based on a mapping 6 : Ca,p(n, (f) — ► W 1 that transforms the cap to 
the unit ball in W 1 according to the following rule: 

S(x) = — — - (x — e n cos <fi) . 
smq) 

Theorem 2. (Levenshtein [T^]). 

rain — 1, — — — - — \ < A(n,8,6) < m(n — 1,2 sin e /2 cot 6). 
V sin q> cos (f> I 

The proof is based on a mapping that projects the cap centrally on the tangent 
hyperplane to the sphere S n ~ x at the point e„. 

Bounds of these two theorems are useful in asymptotics (both as n — > oo and as 
# — > Q) for estimating the size of spherical codes (TO] and the packing density in K™ 
[Ml IT2"] . Their use for finite n is based on the obvious inequality m(n — 1, 2d) < 
A(n, 2 arcsind) and leads to the estimates 

(1) A(n, 8, <p) < A(n + 1,2 arcsin(sin 0/2 cosec <p)) = A(n + 1, lj(6, 4>)) 

(2) A{n,e,4>) < A(n, 2 arcsin(sine/2 cot 0)). 

For large n upper bound ([1]) is uniformly better than bound (J2]) because Ysin^ > 
cot0. For finite n bound ([2]) is stronger than {1} for small <f> and is weaker than it 
otherwise. 



3-B. Spherical strip codes. In this subsection we discuss the action of the pro- 
jection II n on the code in a spherical strip Z(n, 4>])- Given a #-code C C Z 
we would like to know what happens to its distance upon applying the mapping 
II n to it. Given two points xi,X2 G Z the main issue is to establish how the dis- 
tance between their images depends on their relative location in the strip. Let 
dist(xi, X2) — and let the angle 7 7^ (f>, < 7 < 90° be defined by the equation 

cos uj{6, (j)) = coslu(8, (j), 7). 

Geometrically, the angle 7 is defined as follows. Consider two points xi,X2 that 
are 8 away from each other and lie on the boundary of the cap Cap(n, (f) (i.e., 
the angle between each of them and e„ is 0). The distance between their images 
under II n equals lu(9,(/)). Consider the point x' 2 that satisfies Ti n (x2) = ILj^) and 
(x2 7 x' 2 ) = cos 8, then 7 = dist(£ 2 , x' 2 ) (see Figured]). Formally, 7 is the angle given 
by 

sin rWcos 2 8 — cos 2 <f>) 

suit = 

' cos 2 6» + cos 2 0(1 - 2 cos (9) 

Under the mapping H n the code in a strip is transformed to a spherical code in 
5"~ 2 . This transformation yields nontrivial results only in the case of 8 > ip — <p 
(otherwise the code C can contain points that project identically on the equator, 
so the distance of the image code is zero). In this case, the distance in the image 
code is minimized for a pair of points on the "lower" boundary of the strip if 7 < tjj 
and for a pair of points one of which is on the lower and the other on the upper 
boundary, otherwise. 
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n„(x 2 ) ryx,) 

Figure 1. Definition of the angle 7: the triangle x±X2x' 2 is isosceles. 

More precisely, we have the following theorem. 

Theorem 3. (Agrell, Vardy and Zeger pQ) Let < ip < (f> < */a and 9 < 2<p. 
(a) Let 9 > <f> — ip. If 7 < tp then 

A(n,9, =A(n-l,uj(9,4,)). 

If 7 > ijj then 

(3) A(n, 9, </>]) < min {A(n, 0, [if,, 7]) + A(n - 1, 0)), A(n - 1, cj(0, 0, 
(6) Let6 <<j)- ip. Then 

A(n, 6, [if,, </>]) < A(n, 9, ty, 7]) + A(n - 1, </>)) (7 > V 1 ) 
If 9 > 2cj} then A(n, 9, [ip, <j>\) = 1. 

By taking ip = this theorem implies the following corollary for spherical cap 
codes. 

Corollary 1. [1J 

!< A(n,9,-f) + A(n-l,u(9,(f>)) if < 9 < 4> 
= A(n- l,w{6,<t>)) if(f><9<2(t> 
= 1 if 9 > 2<j). 

3-C. Lower bounds. A general lower (existence) bound on the size of a #-code 
in the cap Cap(n, <f>) can be obtained by the standard greedy argument. It follows 
that there exist codes of size 

m > n n ((b)/n n (9), 

where by f2 n (/3) = Y{(n-i)/2) fa sm ™~ 2 T ^ T we denote the area of the spherical cap 
Cap(n,/3) on the sphere S" i_1 . 

4. An extension of Delsarte's method 

In this section we explain a way to use bounds on spherical cap codes in order 
to extend the well-known Delsarte method for bounding the size of spherical codes. 

The original Delsarte (linear programming) bound as applied to spherical codes 
P,[10] has the following form. Let {G^ (x)}k=a,i,... denote the family of Gegenbauer 
polynomials, i.e., polynomials orthogonal on [—1, 1] with weight (1 — x)^" -3 ^ 2 and 
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satisfying the normalization condition G k (1) = 1. Suppose that a real function / 
is a nonnegative linear combination of Gegenbauer polynomials G k n \t), i.e., 

f(t)^Y,^G ( k n \t), where f k > 0. 
fc 

If /0) < for all t £ [-l,cosl?] and f > 0, then A(n,0) < /(l)//o- 

Next we consider an extension of this method to spherical caps. Let Y — 
{yi,...,y m } be a (9-code in the spherical cap Cap(n, cf>) with center e„ and let 
y be the set of all such codes. Of course, m < A(n, 6, </>). Let e* = — e„, let /(£) be 
a real function on the interval [—1, 1], 

H f (Y) = H f ( yi , y m ) := /(l) + f((e*, Vl )) + ... + /«e*, y m », 
h m (n,9,(f),f) ;=m&x{H f (Y)}, h max (n,6,cj), f) := max {h m (n, 6, <j>, /)}. 

m<A(n,8,<f>) 

Theorem 4. Suppose that f is a nonnegative linear combination of Gegenbauer 
polynomials G k (t), i.e., 

/(*) = ^/fc4 n) (*)> where h>0- 
fc 

If f{t) < /or all t £ [— cos 0, cos 0] and /o > 0, then 

At Q\ s h ma,x(n,9,4>,f) 

A(n,0) < . 

Jo 

Proof. Let C — {x\, . . . ,xm} be a 6>-code in S n ~ 1 . It is well known [HI [TO] that 

M M 

^2^2G k n) (ti^) > 0, U,j ■= (xi,Xj) = cos(dist(a;i,x i )). 

i=l 3=1 

Using this we obtain 

M M 

Let 

■7(0 := {J : /«»i,3!i» > 0, j ? »}, = {as, e C : j € J«},m t = 

Note that j G J(i) only if Xj belongs to the Cap(n, 4>) with the center at — x*. 
Then 

Af 

S i (C):=^/« a!i , a;j ))</(l)+ ^ f({x i ,x j )) = Ht{C(i))<h rnt (n,9,<t>,f). 

J=l j6J(i) 

Therefore, 

M 

/oM 2 < S>(C) = ^5i(C) < MVax, 
i=l 

i.e. M < h max / fo as required. □ 

Note that h = /(l). If /(t) < for all t e [-l,cos0], then all m ; = 0, i.e. 
hmax = ho = /(l) and M < /(l)//o, so this theorem includes the Delsarte bound 
as a particular case. 
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For given n, 9, (j), f and m the value h m (n, 9, 4>, f) is the solution of the following 
optimization problem on S n ~ 1 : 

h m (n,9,</>,f) =/(l) + /(- cos <£i) + ... + /(- cos (f) m ) -> max 

subject to the constraints 

<f>i := dist(e* , yi) < <j), 1 < i < m; dist(y i; yj) >9,i^ j. 

The dimension of this problem is m(n — 1) < in — l)A(n,9,(j)). For relatively 
small n and A(n, 9, 4>) optimization can be carried out numerically. Moreover, if 
in addition to the above restrictions the function f(t) is monotone decreasing for 
t S [—1, — cos0] then in some cases the dimension of this problem can be reduced 
to n (see the details in [T71 E]). Suitable polynomials / can be found by linear 
programming (see an algorithm in the Appendix to [14] ^ . 

5. The case of large angles 

In this section we consider 0-codes in a spherical cap Cap(n, (j)) with large values 
of 9. More precisely let us assume that 9 > (p. Clearly, if 9 > 2cf>, then no more 
than one point can lie in Cap(n, </>), i.e. A(n,9,<j)) = 1. Now we consider the case 
2(f>>9>4>. Recall that 4> < tt/2. 

Lemma 1. Suppose 2<p > 9 > <fi > 0, then oj(9,(j>) > tt/3. 
Proof. Let z := cos 9, t := cos cf>. Then 

cos 9 — cos 2 4> z — t 2 z — z 2 z 1 
COSLU{9,(j>) — = < = < 



sm 



\-t 2 ~ \ -z 2 1 + z 2' 



Thus, u(6,4>) > it/3. □ 

As stated in Theorem [3l for 2(f) > 9 > <fi the problem of finding A(n, 9, <f>) is 
equivalent to bounding the size of spherical codes. Since the proof in [1] is not 
isolated into a separate argument we include it here for completeness. 

Theorem 5. (Agrell et al. [J, Musin [13]) If 2c/) > 9 > <j), then 

A(n,9,ct>) = A(n-l,u(6,(t))). 

Proof. First let us prove the lower bound. 

Lemma 2. Let < tp < 4> < tt/2, 9 < 2(f>. Then 

A(n, 9, <t>) > A(n, 9, [ip, <j>\) > A(n - 1, uo(9, 0)). 

Proof. The first inequality is obvious. To prove the second one let consider the 
strip Z(n, [i/j,4>]) C 5 n_1 and let S be its "lower" boundary. The projection n„ is 
a one-to-one map from E to the unit sphere S n ~ 2 (the equatorial sphere of 
Now consider a code C C S n ~ 2 and the code C C E that corresponds to C under 
this map. If the distance of C equals ui(9,4>), then the distance of C is 9 (the 
function uj(9,<fi) is monotone). Since \C'\ = \C\ < A(n, 9, [ip, (/>]), this proves the 
needed inequality. □ 

Now let C = {x±, . . . , x m } be a 9-code in Cap(n, 0). Then 

9i t j :— dist(xi, Xj) > 9 for i ^= j. 

Denote by 4n the angular distance between e n and Xi, where e n is the center of 
Cap(n, </)). Note that <j)i < <fi. 
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Let X — n„(C) be the image of C under the projection on the equator S n ~ 2 
of the sphere from its North pole e„. Denote by 7y = dist(n n (a:,), H n (xj)) be the 
distance between the images of Xj and Xj under the projection. Recall the law of 
cosines for a spherical triangle. Suppose the two sides are a, b and the angle between 
them is ip, then the third side c satisfies 

(5) cos c = cos a cos b + cos tp sin a sin b. 

From this and the inequality cos#jj < cos#, we get 

cos 6i j — cos </> 2 ; cos </>,- cos — cos Si cos </>,- 

cos ^ = • 1 ■ 1 - • — 1 — • — 1 

sm 0i sin (pj sm 0, sin 0j 

^/ cos — cos a cos /? , dQ(a,{3) cos /3 — cos cos a 

Let Q(a,p) = : — , then — = — — . 

sm a sm p da sin a sm p 

From this it follows that if < a, f3 < <b then cos/3 > cos 9 (because #></>); 
therefore dQ(a, f3)/da > 0, i.e., Q(a,f3) is a monotone increasing function in a. 
We have Q(a,/3) < Q(<t>,p) = QifiA) < Q{<j>,$)- Therefore, 

cos 9 — cos d>i cos 6j cos — cos 2 <b 

COS7i,, < ; — - — ; — < r^— = COSUj(9,((>). 

sm 0, sm (pj sm </> 

Thus X is an ui(9, 0)-code on the (n — 2)-sphere. That yields 

A(n,9^)<A(n~l,u(9,(b)). 

□ 

It is proved in |15[ I17j that in the case covered by this theorem, points in an 
extremal configuration are vertices of a convex polyhedron, and lie on the boundary 
of the cap. This implies that if 9 = (f> then the code can be augmented by the point 
e„ without reducing its distance, so A(n,9,9) = A(n,u>(9,0)) + 1. 

Denote by <p n (M) the largest angular distance in a spherical code on S n ~ 1 that 
contains M points. Recall that k(n) denotes the kissing number in n dimensions. 

Corollary 2. Suppose that 9 > 4>, then 

A(n,6,4>) < fc(n-l). 

Moreover, if tp n -i(K) < 7r/3, then 

A(n, 9, 4>) < K. 

Proof. By Lemma Q] we can write uj(9,<j>) = tt/3 + e, e > 0. Then the theorem 
yields 

A(n, 9, (j>) = A{n - 1, u(9, 4>)) = A(n - 1, tt/3 + e) < A(n - 1, tt/3) = k(n - 1). 
If <p n -i(K) < tt/3, then A(n- l,7r/3 + e) < K. □ 

Using Thcorcm[5]together with this corollary we can find the exact value A(3, 9, <fi) 
and A(4, 9, (j>) for > <f>. We can also find A(n, 9, <fi) for all n if cos 9 < cos 2 4>. 

1. Let ?i = 3. Note that fe(2) = 6, ^ 2 (6) = tt/3, and </? 2 (M) = 2tt/M. Then 
A(3,M) = LV^(^^)J <5. 
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2. Let n = 4. In three dimensions the best codes and the values 993 (M) are 
presently known for M < 12 and M = 24 (see [Hlini [22])- It follows from Fejes 
Toth's bound [9 that 

923(2) = 180°, 993(3) = 120°, 993(4) = arccos(-l/3) w 109.47°, 

993(6) = 90°, 993(12) = arccos (1/V5) « 63.435°. 
Schlitte and van der Waerden [22] proved that 
993(5) = 993(6) = 90°, 993(7) « 77.87° (cos 993(7) = cot 40° cot 80°), 

993(8) = arccos — %■ — - « 74.86°, 923(9) = arccos ^ « 70.53°. 
The cases M = 10, 11 were considered by Danzer [BJ: 
993(10) « 66.15°, 993(11) = 953(12). 
Since fc(3) = 12 [23], we have A(4, 9, (/)) < 12. Thus 
A(4, 0, 0) = max {M : 99 3 (M) > w(0, 0)}. 

3. Let cos# < cos 2 </>. In this case we have oj(9,(f>) > 90°. It is well known [21] 
that for all dimensions 

<p n (M) = arccos (jf^j , 2 < M < n + 1; 

and 

9J„(n + 2) = . .. = 9? n (2n - 1) = 99„(2n) = 90°. 
Therefore, for arccos(cos 2 cf>) < 9 < 2(f), we obtain: 

A(n, 9, (/)) — max {M : arccos — > u(6, (f>)}. 

M <n M — 1 

The results obtained can be applied for the kissing number problem as follows. 
For n = 3 let us consider the following polynomial /: 

. 2431 q 1287 7 18333 » 343 , 83 , 213 9 f 1 

f(t) = 1 9 1 7 H 1 5 H 1 3 i 2 H . 

3 w 80 20 400 40 10 100 10 200 

This polynomial satisfies the assumptions of Theorem 2] with 9 — 60°, (j> « 53.794° 
(/(— cos0) = 0), and /o = 1. In this case A[$,0,<j>) = 4. Since 7i max < 13 (see a 
proof in [IS]) we have jfc(3) = A(3,7r/3) < 13, i.e. fc(3) = 12. 
In the case of n = 4, Theorem [4] can be applied with 

/(f) = 53.76i 9 - 107.52i 7 + 70.56t 5 + 16.384i 4 - 9.832t 3 - 4.128t 2 - 0.434< - 0.016. 

Here 9 = 60°, <j> « 52.559°, / = 1, and A(4,9,<p) = 6. It was proved [H] that 
ft max < 25. Since fc(4) > 24 this yields fe(4) = 24. 

Recently, Pfender [20] considered the case cos 9 < cos 2 (f>. He found some im- 
provements for upper bounds on k(n) for dimensions n = 9, 10, 16, 17, 25, 26. 
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6. Stretching transformation 

In this section we will prove the following bound on spherical cap codes which 
relates the maximum size of such a code to the size of codes in a hemisphere. 

Theorem 6. Let e / 2 < (ft < V 2 - Then 

(6) A(n,e,</>)<B(n,w(0,</>)). 

Remark. This theorem improves upon the bound (fTJ) by reducing the dimension 
on the right-hand side by one. It also extends the applicability of the bound in 
Theorem [5] to the range of angles 8/2 < (ft < 8, although in this range we cannot 
claim the exact equality anymore. On the other hand, by Theorem [6] it is sufficient 
to estimate the number of code points in the hemisphere as opposed to the entire 
sphere. 

Note also that for any 8, (ft, 

cosui(8, (ft) — cos 8 — tan 2 <ft(cos8 — 1) < 

therefore, this theorem is stronger in the entire range of angles than the trivial 
bound A(n, 8, (ft) < A(n, 8). Finally, the angle uj(8, (ft) ranges between 7r and 8 as (ft 
grows from e /2 to w /2 and is a monotone decreasing function of 8. 

A proof of Theorem [6] will follow from the following result which describes the 
effect on the distance of spherical cap codes of a "stretching map" of spherical caps. 

Theorem 7. Let < (ft < n/2. Then for any s > 1 

A(n,e,(j>) < A(n,8',s(ft), 

where 

ai 2 , sm2 s / n 2 i \ 
cos 8 — cos scp H 5 — (cos 8 — cos <p). 

sin (j> 

Proof: Let T s ,s > 1 be a map on Cap(n, <j)) defined as follows: for a point 
x G Cap(n, (f>) with dist(cc, e n ) = a its image y = T s (x) is a point that satisfies 
dist(y, e n ) = sa and lies on the meridian passing through x. Thus, T s (Cap(n, </>)) = 
Cap(n, s4>), and we assume that s<p < V 2 - 

The proof of Theorem [7] relies upon the next two lemmas. 

Lemma 3. Let x\,X2 G Cap(n, (ft) with dist(xi, e„) — u,dist(x2,e n ) = v, 
dist(xi, X2) = 8. The distance dist(T s (xi), T s {x2)) reaches its minimum when u = v. 

Proof. Figure [2] shows the relative location on the sphere of x\, a; 2 and their images 
yi = T s (xi),y 2 = T s {x 2 ). 




Figure 2. 
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Using ([5]) we find that cos 9' — F(u, v, 8), where 

F(u,v,6) — cossucossv + p(u)p(v)(cos9 — cos u cos v), 

p(t) = sin(si)/ sin t. For definiteness assume that v < u. We need to prove that 
F(u, v, 9) < F(u, u, 9) where < u — v < 9 and su < tt/2. 

Fact 1. (i) The function p(t) is monotone decreasing for t 6 (0,7r/2s). Indeed 

, cos t cos st 

P V-) = 2 \ s ^ an * — tan st) < 

sin t 

with the equality only if t = 0. 

(ii) /?(£) < s (follows from (i) and the equality p(0) = s). 

Fact 2. The function S(9) = w, 9)—F(u, u, 9) is maximized on 9 for = u— v. 
Proof: The coefficient of cos6> in S(9) equals p(u)(p(v) — p(u)) > 0. Then the 
claim follows from the condition ^ji > 9 > u — v. 

Fact 3 (which implies the lemma). S(9) < S(u — v) < 0. 
Proof: The first inequality is proved in Fact 2. Now compute 

■ 2 

S(u — v) = cos s(u — v) — cos 2 su ^ — (cos(m — v) — cos 2 u) 

sin u 

= cos s(u — u) — 1 — p 2 (u) cos(m — v) + p 2 {u). 
The derivative of the last expression on v equals (sp(u — v) — p 2 (u)) sin(u — v). Since 
p{ u ) 5: p{ u ~ v), we can write 

{sp(u — v) — p 2 {u)) sin(w — w) > (s — p{u))p{u — v) sin(u — v) > 
where the last inequality follows from part (ii) of Fact 1 above. □ 

Lemma 4. Letx\,X2 € Cap(n, </>) wii/i dist(xi, e n ) = dist(a;2, e„) = u, dist(xi, x-i) 
II . where u < <j>. Then the distance dist(T s (xi), T s (x2)) reaches its minimum when 
u = <f>. 

Proof. Since 

cos6»' = F(u,u,9) = 1 - p 2 (u){l - cos0), 
the claim is implied by Fact l(i) above. □ 

The last two lemmas imply Theorem [7] Indeed, let x\,X2 € Cap(n,#, </)) be two 
points at distance 9. The lemmas show that in order for the distance between their 
images under T s to reach its minimum the points should lie on the boundary of the 
cap. Then the expression for 9' in the theorem is implied by an application of the 
cosine law J5]). 

Finally, Theorem [6] follows from Theorem [7] by taking s — 7r/2</>. □ 
We conclude this section with two applications of Theorem [6] 

6- A. An upper bound on spherical codes. Here we establish the following 
new estimates: 

(7) A {n,6) < ^-B{n,u{0A))- 

(8) A (n,6)<-^-A(n-l,u;(9,4,)) (9 > <j>) 

where f2 n (</>) is the area of the cap of radius <f> and Q n = n n / 2 /T( n /2 + 1) is the 
"surface area" of the unit sphere S'' 1 ^ 1 . They are implied by the Bassalygo-Elias 
inequality stated in the next lemma. 
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Lemma 5. HE] Let e /2<4>< */2. Then 

A{n,0) < J*»A(n,9,<l>)' 
u n {(p) 

Proof: Consider a code C C S" 1-1 and let C<p(z) be the number of code points 
in the cap with "center" z and radius (f>. Note that every cap whose center z is at 
most 9 away from a given code point x will contain this point. Then clearly 

/ \c${z)\dz = o n (0)|q. 

Since C^,(z) < A(n,9,<fi), we obtain 

A(n,9,^)a n > A(n,9)n n {q>). 

□ 

Therefore, using Theorems [5] and [5] we obtain the bounds J7|), ©. In particular, 
inequality (JT]) is stronger than bounds on cap codes based on Lemma[5]that appeared 

in [2i que]. 

6-B. Large dimensions. Let R{C) = l /n\n \C\ be the rate of the code C C S n ^ 1 . 
Denote by 

R + {6) =limsup-lnA(n,0) R~{9)= liminf - ]nA(n, 9). 

n-i-oo n n->oc n 

Abusing notation, we write below R(9) to refer to the common value of R + and R~ 
even though it is not known that the limit exists. Likewise we write R{9, <f) and 
R(9, [V>, 4>]) to refer to cap and strip codes. In this section we show that the problem 
of finding either of the last two quantities is equivalent to that of computing R{9). 
More precisely, we have the following theorem. 

Theorem 8. Let < ip < (f> < tt/2. Then 

R(9, hM) = R{9,4>)=R{u{9,4>))- 

Indeed, combining Lemma [21 Theorem [5] and the obvious inequality B(n, a) < 
A(n, a) we obtain the bounds 

(9) A(n - 1, u{9, <{>)) < A{n, 9, [>ip, 0]) < A{n, 9, 0) < A{n, lo(9, </>)) 

which imply this theorem. 

Inequality ^ can be rewritten in a somewhat more visual way. Consider two 
points X\, X2 G Cap(n, <p) with angular distance 9. Let d be the Euclidean distance 
between them. From ([5]) the angular distance between their images under II n is 9' < 
u>(9,(f>) and the Euclidean distance equals 2sin e '/2 < d/ sin0. Hence the minimum 
distance of the image code II n (C) c S n ~ 2 is at most dj sin^. Denote by N(n, d, (f>) = 
A(n, arccos(l — d /2),<p) the maximum number of points in a spherical cap code 
with minimum Euclidean distance d and let N(n, d) be the same for the sphere. An 
equivalent form of ([9]) for spherical caps is as follows: 

N(n- l,d/sin0) < N(n,d,<f>) < N(n,d/ sm<j>). 

A somewhat weaker upper bound on A(n, 9, 4>) than in ^ is given by fl}. This 
bound is nevertheless sufficient to establish one part of the asymptotic claim of 
Theorem HI 
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7. Codes in hemispheres 

Here we consider upper bounds on B(n,9) = A[n,9 ,tt /2). Let C C 5+ = 
Cap(n, 7r/2) be a #-code in the hemisphere. Denote 

(10) C{[a,0\) =CnZ(n,[a,/3]), C{a) = C n Cap(n, a). 

Theorem 9. Let 9 < */2, 5 = (n — 6)/2 and let C C S + be a 9-code. Then 

\C{[6,-K/2])\ + 2\C{5)\<A{n,e). 

Proof. Let a — \C([S, 7r/2])|, b = \C(S)\. For a point x — (x%, . . . , x n -%, %n) denote 
by x* — (xi, . . . , x n -Xi — x n ) its reflection about the equator. Let C*(5) — {x* : 
x e C(S)} be the reflection of the code C(5). Consider the code Q — C U C*(5). 
We claim that Q is a 6>-code. Referring to Fig. [3] this amounts to showing that 
dist(q,p*) > 6 if dist(q,p) > 0. To prove this, we choose the point s so that the 
angle Zqsp = 90° and use (|5|) as follows: 

cos 9 > cos a = cos rj cos /3 > cos r] cos /?* = cos a* . 

This proves that Q is a 6»-code. Then a + 2b = \Q\ < A(n, 9). □ 

Corollary 3. 

B(n, 9) < i/a(A(n, 0, [(vr - 6»)/2, tt/2]) + A(n, 9)) 
Proof. Using the notation of the previous theorem, we have 

2\C\ = 2a + 2b<a + A{n, 9) < A(n, 9, [S, tt/2]) + A(n, 9). 

□ 




Figure 3. 

Corollary 4. 

^ ^ A(n-1,6)+A(n,e) - cos 9 

(11) B(n,9)< , cos0 = — _. 

Proof. We use Theorem [3] In our case <f> — tt/2, ip = (it — 9)/2, therefore, 9 > 
<j> — i/j > 9/2, so part (a) of this theorem applies. Substituting the values of <j>,ip 
in the inequality A(n, 9, [ip, <j)]) < A(n — 1, uj(9, <p, tjj)), we obtain A(n, 9, [i/j, tt/2]) < 
A(n — 1,9). The proof is concluded by using this estimate in Corollary [3] □ 

Using the value 9 — 60°, we obtain estimates on the one-sided kissing number 
B(n). In particular, since A(n, 60°, 30°) = 2, the last corollary yields 

Advances in Mathematics of Communications Volume 1, No. 1 (2007), 1-19 



14 



A. Barg and O. R. Musin 



Corollary 5. 

B(n) < min[(l/2)(A(n- 1,770) +k(ri)), k(n) — 2], rj := arccos w 54.74°. 

V3 

Clearly, 5(2) = 4. Let us use this bound for n — 3,4. 

n = 3. We have fc(3) = 12, >i(2, 770) = 6, then 5(3) < 9. On the other hand, 
5(3) > 9, so -B(3) = 9. Note that in this case the bound is sharp. 

n = 4. Recently, K. Bezdek [21 [3] proved that 5(4) — 18 or 19, and conjectured 
that 5(4) = 18. It was proved, also recently [HUH], that fc(4) = 24. Delsarte's 
linear programming method gives A(3, 770) < 15. Thus 5(4) < 19. The proof that 
5(4) = 18 given in [T5] is based on an extension of Delsarte's method. 

For higher dimensions we can rely on the known bounds for spherical codes. 
Denote by g n an upper bound on A(n,7?o) given by Delsarte's linear programming 
method, and by k n the known upper bounds on k(n) (see, e.g., Table 1.5 in [S]). 
Then Corollary [5] implies 

tdi \ ^ 9n-i + k n 
B(n) < . 

This gives the following bounds: 

5(5) < 39, 5(6) < 75, 5(7) < 135, 5(8) < 238. 

For instance, for n = 8 we have gj = 236 (obtained by Delsarte's method with a 
polynomial of degree 11) and k$ = 240. Note that g n -i > kn for n > 8; therefore, 
for these dimensions we just have the bound B(n) < k(n) — 2. 

However, even for 5 < n < 8 these bounds are not sharp. Our conjectures for 
n = 5, 8 are 

5(5) = 32, 5(8) = 183. 

8. A BOUND ON SPHERICAL CAP CODES 

In this section the methods and results developed above will be used to derive 
another bound on spherical cap codes. Given a code C C Cap(n, 0), our plan is to 
first map the cap on the hemisphere S+ and then use the results of the previous 
section together with some additional ideas. In the next theorem we use notation 

do). 

Theorem 10. Let < <fi < w/2,u)x := u>(6,<p),u)2 '■= <^{0,<p,tp), where ip — 0(1 — 
W1/71"). Then for any 9-code C C Cap(n, 4>) we have 

(12) |C(hM)| < A(n-l,min(o;i,W2)) 
and 

(13) |C([^])|+2|Cty)|<A(n,wi). 

Proof. Consider the mapping T s ,s — n/2(j) that sends the cap to the hemisphere 
S+. By Theorem [71 the code C is mapped to a code T S (C) C 5+ with distance io\. 
Now (fT3]) is implied by Theorem [9) 

To prove (fl"2"]) . let us bound above |C([^),0])| in terms of 9 and <fi. Consider 
the action on the code C([ip,4>]) of the orthogonal projection LI„ on the sphere 
from e n on the equator. The code D = U n (C([tp,(j)])) is a spherical code in 
n — 1 dimensions. Given two points x,y € C([/0, </>]) such that dist(x,e n ) = 
«i, dist(i/, e„) = ot2, dist(x, y) — (3, the distance between their images under II n 
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is given by a\, 0.2 ) ■ As shown in the proof of Lemma[2j this function is mono- 
tone on ot\ (or 0L2 ) if the other two arguments are fixed, so its minimum is attained 
at one of the boundaries. Therefore, the distance of the code D is determined 
according to one of the following two cases: 

(i) There are two points C\ , C2 G C lying on the boundary of the cap and 9 away 
from each other. Their images in the code D are two points on the equator 
at distance u>\. 

(ii) There are two points ci, C2 in C at distance 9 such that T s (c\) — x%, T s (c2) = 
x 2 (see Fig. |4j). 















/v \ 






/ c 






31-03, 

2 / 




/ 


' c \ 
1 \ 









Figure 4. 

Upon projecting X2 on the equator, the distance dist(II n (x2), x\) = u>(6, <f>, ip), 
where ip = dist(e„, C2) satisfies 

7T - Wl 

i.e., dist(n„(x 2 ),xi) = lu 2 . 
Therefore, D is a code with distance at least min(iJi, W2). □ 

Corollary 6. For all < 9 < 4> < V 2 

(14) A(n, 0, 0) < l /2(A(n - 1, min(wi, w 2 )) + A(n, wi)). 

Proof. Use (| 12[) and (fl~3| to compute 

2|C(hM])|+2|<7(0)| = |C([^,^])|+|C(^)|+|C(^)| <A(n-l,miix(w 1 ,W2))+A(n,«i). 
This gives the claimed result. □ 

Comparison of this bound with the other bounds considered in this paper is 
difficult in general. However, it is clear, that it improves upon the bound ([6|) in all 
cases when ui\ < lo2- The domain of values of 9, </3 for which this holds true is shown 
in Fig. [5] 
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Figure 5. Area of improvement for bound (|14jl 

Examples. Next we give several examples of the bounds known earlier and 
obtained in this paper, showing that the results of this paper improve the known 
estimate in some range of values of 9 and <f>. 

Bounds on A(6, 9, 4>) 



9 


4> 


Corollary Q] 


Theorem [H] 


Corollary [5] 


0.027T 


0.27T 


9069268 


7581886 


4058040 


0.057T 


0.25tt 


145587 


194908 


111241 


0.2vr 


0.47T 


661 


832 


591 


0.25tt 


0.47T 


221 


272 


201 


0.37T 


0.47tt 


174 


138 


115 



Bounds on log 2 A(n, 9, 4>) 



n 


9 





Corollary [T] 


Theorem [6] 


Corollary [6] 


20 


0.057T 


0.47T 


66.1912 


65.7786 


65.3665 


20 


0.087T 


0.47T 


52.7212 


52.891 


52.6557 


30 


0.057T 


0.47T 


97.75 


97.7008 


97.8523 


40 


0.057T 


0.45tt 


132.516 


131.447 


133.45 


100 


0.057T 


0.45tt 


322.205 


321.917 


333.28 



To compute the bounds in examples one needs to use some upper bound on 
A(n,9). We have used the bound of Levenshtein [12], [13l p. 618] which is the best 
known universal upper bound on spherical codes. 

9. A GENERAL METHOD OF BOUNDING THE SIZE OF CAP CODES 

In this section we generalize the method of Sections [JJ [5] to develop a general 
approach to bounding the size of codes in spherical caps. 

Let C £ Cap(#, 0) be a #-code. Let = 0o < 4>i < • ■ ■ < 0fc < ^fc+i = <f>- 
To bound above the size of the code C we partition the cap as follows: 

Cap(n, 4>) = Cap(n, <j> x ) U Z(n, [<t>\,<fa]) U • • • U Z(n, [ipkAk+i])- 
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Let pi := \C n Z(n, [^i>^i+i])l be the size of the code in the strip which can be 
estimated by the methods of [T] and this paper. In particular, denote by a$j an 
upper bound on the size of the code in the strip Z(n, [4>i, <fij]) and let a, := o^t+i. 
Suppose that we can compute these estimates for some subset P of pairs We 
have the following linear constraints. 

(15) 0<Pi<di, (i,i + l)eP; Pi + . . .pj-i < dij, (i,j)eP. 

Another set of linear inequalities can be obtained from the arguments similar to 
the proof of Theorems [9l [TOl Namely let r(i) — max{j : s<fij < (tt — ui(6, <j>i+i))/2}. 
Consider the codes C = C(</V(;)+i) an d C" — C([0 r (j) + i, (pi+i])- Upon stretching 
the cap to the hemisphere by applying the mapping T s and using (13)) , we obtain 
2|C"| + \C"\ < bi, where hi is the size of spherical code with distance ui(6, (f> r ^ + i). 
We then obtain inequalities 

(16) 2(p Q +p 1 -\ \-Pr(i)) +Pr(i) + l H \-Pi<h, 

where $ is a subset of indices for which we perform the described procedure. 
We summarize the arguments of this section in the following 

Theorem 11. The size of the code C is bounded above by the solution of the linear 
programming problem 

Po+Pi^ h Pk —> max 

under the constraints H5\) . U6\) . 

Note that Theorem [10] constitutes a solution of this problem in the case k = 0. 

The set {4>i} can be optimized in computations. It can be formed, for instance, 
by taking all the angles 7 produced by the recurrent calculation in Corollary [1] 
together with some additional breakpoint angles. 

10. Applications 

We have discussed the applications of the bounds on codes in caps to the kissing 
number and the one-sided kissing number problems. As remarked in the introduc- 
tion, the methods developed in this paper can be also useful in the problems of 
estimating the packing density in M. n and of the size of constant weight codes. We 
end the paper with brief remarks on these applications. 

10-A. Spherical codes and packing density. Let A„ be the density of packing 
the n-dimensional real space with equal nonoverlapping balls. A classical problem 
in geometry is to compute A„ for a given n and for n — > 00. It is known 5, p. 265] 
that 

(17) A„ < i/2(sin s/ 2 ) n A(n + 1,6) (0 < 6 < tt). 
Using Lemma [5] we now obtain the estimate 

(18) A w < 1/2 (sin °/2) n Q n+1 A( " + ^ (0<9< n,6/2 < < tt/2). 

To compute upper bounds on A„ for a given value of n we have a choice of using (TT7)) 
with the known bounds on A(n + 1, 0), or using (fT8)) together with the bounds on 
spherical cap codes considered in this paper. Observe that the best known estimate 
of the packing density for n — > 00 [10] is obtained by employing (fT8| together with 
the bound (fT]) on cap codes. Therefore, inequality ([18]) coupled with better bounds 
on cap codes derived in this paper will also improve the density estimates for finite 
(but possibly large) values of n. 
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10-B. Constant weight codes. A constant weight binary code is a subset of 
{0, 1}™ formed of vectors with a fixed number, say w, of ones. Denote by A(n, d, w) 
the maximum size of a constant weight code with minimum Hamming distance d. 
Computing or estimating the numbers A(n, d, w) is a problem with a long history 
in coding theory, summarized in pQ. The most studied region of parameters is 
n < 28 for which the most recent tables are published in pQ. For larger n bounds 
on A(n, d, w) are tabulated in |25j . In [1] the problem of bounding above A(n, d, w) 
was reduced to bounds on spherical cap codes which led to several improvements of 
the tables for short lengths. We intend to use the new bounds on cap codes derived 
in this paper to further improve the tables. 

As a final remark, note that all the bounds on cap codes considered in this paper 
rely on bounds on the maximum size A(n,0) of spherical codes. In calculations, 
apart from the Levenshtein bound and related results it is possible to use a direct 
solution of Delsarte's linear programming problem relying on a method developed 
in [19]. 

<abarg@umd. edu, omusinOgmail . com> 
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